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Abstract: We consider second-order viscous hydrodynamics in conformal field theories at 
finite temperature. We show that conformal invariance imposes powerful constraints on the 
form of the second-order corrections. By matching to the AdS/CFT calculations of correla- 
tors, and to recent results for Bjorken flow obtained by Heller and Janik, we find three (out 
of five) second-order transport coefficients in the strongly coupled N = 4 supersymmetric 
Yang-Mills theory. We also discuss how these new coefficents can arise within the kinetic 
theory of weakly coupled conformal plasmas. We point out that the Miiller-Israel-Stewart 
theory, often used in numerical simulations, does not contain all allowed second-order terms 
and, frequently, terms required by conformal invariance. 
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1. Introduction 



Relativistic hydrodynamics is an important theoretical tool in heavy-ion physics, astrophysics, 
and cosmology. Hydrodynamics gives reliable description of the non-equilibrium real-time 
macroscopic evolution of a given system. It is an effective description in terms of a few 
relevant variables (fields) and it applies to the evolution which is slow, both in space and in 
time, relative to a certain microscopic scale ]l[ |2|. 

In the most common applications of hydrodynamics the underlying microscopic theory 
is a kinetic theory. In this case the microscopic scale which limits the validity of hydrody- 
namics is the mean free path £ m f p . In other words, the parameter controlling the precision of 
hydrodynamic approximation is k£ m f p , where k is the characteristic momentum scale of the 
process under consideration. 

More generally, the underlying microscopic description is a quantum field theory, which 
might not necessarily admit a kinetic description. An experimental example of such a system 
is the strongly coupled quark-gluon plasma (sQGP) recently discovered at the Relativistic 
Heavy-Ion Collider (RHIC) at Brookhaven National Laboratory. The N = 4 supersymmetric 
SU(N C ) Yang-Mills theory in the limit of strong coupling provides a theoretical example of 
such a system which, in the limit of large number of colors N c , can be studied analytically 
using the AdS/CFT correspondence ||. In these cases, where kinetic description may be 
absent, the role of the parameter £ m f p is played by some typical microscopic scale. In the 
above examples this scale is set by the temperature: £ m { p ~ T . 

When the parameter k£ m { p is not too small, one may want to go beyond the first order 
in k£ m f p . This is the case, for example, in the early stages of heavy-ion collisions. There are 
two sources of corrections beyond the k£ m f p order. First, there are corrections due to thermal 
fluctuations of hydrodynamic variables contributing via nonlinearities of the hydrodynamic 
equations. The fluctuation corrections lead to nonanalytic low-momentum behavior of certain 
correlators [Q] (similarly to the chiral logarithms that emerge from loops in chiral perturbation 
theory) and are, for example, essential for describing non-trivial dynamical critical behavior 
near phase transitions || . Such corrections are calculable in the framework of hydrodynamics 
with thermal noise. 

The second source of corrections are second-order terms (order (k£ m f p ) 2 ) in the hydrody- 
namic equations, sometimes called the Burnett corrections ||. These corrections come with 
additional transport coefficients. These second-order transport coefficients are not calculable 
from hydrodynamics, but have to be determined from underlying microscopic description or 
input phenomenologically, similarly to first-order transport coefficients such as viscosity. 

In gauge theories with a large number of colors N c the corrections of the first type 
(fluctuation) are suppressed by 1/N% an< ^ therefore the corrections of the second type 
(Burnett) dominate in the limit of fixed k and N c — ► oo. For this reason, in this paper, 
we concentrate on the second type of corrections. Moreover, we shall consider the case of 
conformal theories, where the number of second-order transport coefficients is substantially 
reduced. In the real-world applications we deal with fluids which are not exactly conformal, 
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however, e.g., QCD at sufficiently high temperatures is approximately conformal. 

This paper is organized as follows. In Sec. § we derive the consequences of conformal 
symmetry for hydrodynamics. In Sec. [| we classify all terms of order k 2 consistent with 
conformal symmetry. In Sec. Q we compute three of the five new transport coefficients for 
the strongly-coupled M = 4 supersymmetric Yang-Mills (SYM) theory using the AdS/CFT 
correspondence. In Sec. [| we show that hydrodynamic equations derived from the kinetic 
description (Boltzmann equation) of a weakly coupled conformal theory do not contain all 
allowed second-order terms. In Sec. |(| we analyze our findings from the point of view of 
the Miiller- Israel-Stewart theory |7|, ||, ||, [10[] , which involves only one new parameter at the 
second order, and show that this parameter cannot account for all second-order corrections. 
Our conclusions are summarized in Sec. 0. 



2. Conformal invariance in hydrodynamics 



To set the stage, let us emphasize again that hydrodynamics is a controlled expansion scheme 
ordered by the power of the parameter k£ m fn, or equivalently, by the number of derivatives 
of the hydrodynamic fields. We shall set up this expansion paying particular attention to the 
consequences of the conformal invariance on the equations of hydrodynamics. 

2.1 Conformal invariance and Weyl anomalies 

The hydrodynamic fields are expectation values of certain quantum fields, such as e.g., com- 
ponents of the stress-energy tensor, averaged over small but macroscopic volumes and time 
intervals. Such averages can, in principle, be calculated in the close-time-path (CTP) formal- 
ism [11]. Consider a generic finite-temperature field theory in the CTP formulation. Turning 



on external metrics on the upper and lower contours, the partition function is 

z \9l*>9%] = j V4>iV4> 2 exp {iS[(j)i, g^] - iS[fc, g^]} , (2.1) 

where cfii and <p2 represent the two sets of all fields living on the upper and lower parts of the 
contours, and S[(f>, g^ u ] is the general coordinate invariant action. 

The one-point Green's function of the stress-energy tensor is obtained by differentiating 
the partition function (the metric signature here is — h ++): 

(T'n = - _ , 1 , 2.2 
(T 2 ^) = -=— -5-, (2.3) 



V-92 $g 2 



where (. . .) denote the mean value under the path integral and y/—gi t 2 = y — detg^. 

In this paper we consider conformally invariant theories. In such theories the action 
S[<t>,gfj,u] evaluated on classical equations of motion 5S/5(j) = and viewed as a functional of 
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the external metric g^ v is invariant under local dilatations, or Weyl transformations: 



9pv -> e ^g^, (2.4) 

with parameter u a function of space-time coordinates. As a consequence, classical stress- 
energy tensor T^ u = SS/Sg^ is traceless since g^T^ = —(l/2)5S/5u> = 0. 
In the conformal quantum theory (^4]) the Weyl anomaly []i~2| , |i~3|| implies 

9UT l n = W d [g] lu l (2.5a) 
9 UT 2 n=W d [g%], (2.5b) 

where W d is the Weyl anomaly in d dimensions, which is identically zero for odd d. For d = 4: 

W 4 [^] = -^(R^xpR^-^uR^+R^ + r^iR^xpR^-^R^R^ + ^R 2 )^ (2.6) 

where R^uXp and i?^ (-R) are the Riemann tensor and Ricci tensor (scalar), and for SU(N C ) 
M = 4 SYM theory a = c=\ (iV c 2 - l) Q. The right-hand side of Eqs. ((27^) contains four 
derivatives. In general, for even d = 2k, W2k contains 2k derivatives of the metric. 

Let us now explore the consequences of Weyl anomalies for hydrodynamics. The hydro- 
dynamic equations (without noise) do not capture the whole set of CTP Green's functions, 
but only the retarded ones. Hydrodynamics determines the stress-energy tensor T^ v (more 
precisely, its slowly varying average over sufficiently long scales) in the presence of an arbi- 
trary (also slowly varying) source g^ v . The connection to the CTP partition function can be 
made explicit by writing 

d\u = 9pv + ]ppv, 9%, = 9pv ~ \ppv (2.7) 

If Ipv = then Z = 1 since the time evolution on the lower contour exactly cancels out the 
time evolution on the upper contour. When 7^ is small one can expand 

lnZ = ^ Jdx^gJx) lllu (x)T^(x)+0( 1 2 ), (2.8) 

where T fiu (x) depends on g^ u , and is the stress-energy tensor in the presence of the source 
gnv At long distance scales it should be the same as computed from hydrodynamics. 

Substituting Eqs. (|^) and (H) into E Q- dU)> the 0(1) and 0(7) terms yield two 
equations: 

9fXU T^ = W d [g^}, (2.9a) 



{x) S[y^)Tpx)] + ^—^ Wd[g (x)]) . 

(2.9b) 

In odd dimensions, the right hand sides of Eqs. (|2.9j ) are zero. In even dimensions, they 
contain d derivatives. In a hydrodynamic theory, where one keeps less than d derivatives, 
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they can be set to zero. For example, at d = 4, the Weyl anomaly is visible in hydrodynamics 
only if one keeps terms to the fourth order in derivatives. This is two orders higher than 
in second-order hydrodynamics considered in this paper. For larger even d, one has to go 
to even higher orders to see the Weyl anomaly. Thus, we can neglect Wd on the right hand 
side: second-order hydrodynamic theory is invariant under Weyl transformations. The two 
conditions ( |2.9D then become 

9lxv T^ = 0, (2.10) 



^w = -[i +1]6{x - y)T {x) - (211) 



Since the r.h.s. of equation ( p. 11 ) is — (1/2)5T^ U /Suj it implies the following tranformation 



law for T^ v under Weyl transformations (|2.4| ): 

rpflU _^ e {d+2)uj rpflU^ (2.12) 

Noting that In Z is invariant under Weyl transformations this could have been gleaned from 
Eq. (^^) already. 

A simple rule of thumb is that for tensors transforming homogeneously 

^m—Hm _^ e A,4u^(ii.../i m (2.13) 

the conformal weight Aa equals the mass dimension plus the difference between the number 
of contravariant and covariant indices: 



Aa = [A]+m- n. (2.14) 
2.2 First order hydrodynamics as derivative expansion 

The existence of hydrodynamic description owes itself to the presence of conserved quantities, 
whose densities can evolve (oscillate or relax to equilibrium) at arbitrarily long times provided 
the fluctuations are of large spatial size. Correspondingly, the expectation values of such 
densities are the hydrodynamic fields. 

In the simplest case we shall consider here, i.e., in a theory without conserved charges, 
there are 4 such hydrodynamic fields: energy density T 00 and 3 components of the momentum 
density T°\ It is common and convenient to use the local velocity instead of the momentum 
density variable. It can be defined as the boost velocity needed to go from the local rest frame, 
where the momentum density T° l vanishes, back to the lab frame. Similarly, it is convenient 
to use e - the energy density in the local rest frame - instead of the T 00 in the lab frame. 
The 4 equations for thus defined variables e and are conservation equations of the energy- 
momentum tensor VuT^" = 0. 

In a covariant form the above definitions of e and can be summarized as 

T^ v = eu^u v + Tj"\ (2.15) 
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In hydrodynamics, the remaining components (spatial in the local rest frame: u^T 1 ^ = 0) 
of the stress-energy tensor T^ u appearing in the conservation equations are not independent 
variables, but rather instantaneous functions of the hydrodynamic variables e and and 
their derivatives. In the hydrodynamic limit, this is the consequence of the fact that the 
hydrodynamic modes are infinitely slower than all other modes, the latter therefore can be 
integrated out. All quantities appearing in hydrodynamic equations are averaged over these 
fast modes, and are functions of the slow varying hydrodynamic variables. The functional 
dependence of (constituitive equations) can be expanded in powers of derivatives of e 
and u^. 

Writing the most general form of this expansion consistent with symmetries gives, up to 
1st order in derivatives, 

Tf = P(e)bT - 7]{e)a^ - ({e)A^{V-u), (2.16) 
where the symmetric, transverse tensor with no derivatives is given by 

In the local rest frame it is the projector on the spatial subspace. The symmetric, transverse 
and traceless tensor of first derivatives a^ v is defined by 

a*" = 2<VV> , (2.18) 

where for a second rank tensor the tensor defined as 

<^> = Ia^A^A^ + A Pa ) - ^-j-A^A Q/3 A a/3 = (2.19) 

is transverse u^A^^ = (i.e., only spatial components in the local rest frame are nonzero) 
and traceless g^ u A^ = 0. 

In the gradient expansion fl2.16D , the scalar function P(e) can be identified as the ther- 
modynamic pressure (in equilibrium, when all the gradients vanish), while 77(e) and £(e) are 
the shear and bulk viscosities. The expansion coefficients P, rj and £ are determined by the 
physics of the fast (non-hydrodynamic, microscopic) modes that have been integrated out. 

2.3 Conformal invariance in first-order hydrodynamics 

It is straightforward to check that if T^ v transforms as in Eq. fliTTD and = 0, then its 



covariant divergence transforms homogeneously: V^T^ — ► e( d+2 ^V M T^, hence the hydro- 



dynamic equation V^T^ = is Weyl invariant [15]. 

Let us now see what restrictions conformal invariance imposes on the first-order consti- 
tutive equations ( p.16 ). First, the tracelessness condition Tff = forces e = (d — 1)P and 



C = 0. Since in a conformal theory e = const-T rf , we shall trade e variable for T in what 
follows. Since g^ u u^u u = — 1 the conformal weight of is 1. By definition ( |2.15| ) and by 
( |2.12j ) e has conformal weight d and therefore 

T -» e^T, u» -> e^u^ (2.20) 
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in accordance with the simple rule ( |2.14| ). 
By direct computation we find that 



_» e^a^, (2.21) 

i.e. cP" transforms homogeneously with conformal weight 3 independent of d (in agreement 
with ( 2.14[ )). For conformal fluids 77 = const • T d ~ l , and therefore T^ v transforms homoge- 



neously under Weyl transformation as in Eq. ( |2.12| ). 
3. Second-order hydrodynamics of a conformal fluid 

In this Section we shall continue the derivative expansion ( |2.16| ). We shall write down all 
possible second-order terms in the stress-energy tensor allowed by Weyl invariance. Then we 
shall compute the coefficients in front of these terms in the M = 4 SYM plasma by matching 
hydrodynamic correlation functions with gravity calculations in Section ||. 

3.1 Second-order terms 

Rewriting Eq. ( 2.15| ) we introduce the dissipative part of the stress-energy tensor, IP": 

TV = eu^v? + PA^ + IP" , (3.1) 

which contains only the derivatives and vanishes in a homogeneous equilibrium state. The 
tensor IP" is symmetric and transverse, u^IP" = 0. For conformal fluids it must be also 
traceless g lw W w = 0. To first order 

IP" = -r?cP" + (2nd order terms), (3.2) 



where cP" is defined in Eq. (2.18|). We will also use the notation for the vorticity 



jy* = Ia^A^V^ - V p u a ) . (3.3) 

We note that in writing down second-order terms in IP", one can always rewrite the 
derivatives along the d- velocity direction 

D = n^V M (3.4) 

(temporal derivative in the local rest frame) in terms of transverse (spatial in the local rest 
frame) derivatives through the zeroth-order equations of motion: 

D1nT = — -!— (Vx-tt), Du" = -Vf InT, Vf = A^ Q V Q . (3.5) 
d — 1 

Notice also that V 1 • u = V • u. 
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With the restriction of transversality and tracelessness, there are eight possible contribu- 
tions to the stress-energy tensor: 



V^lnTVhnT, V^V^lnT, a^(V-u), a^ x a 



a 



(3.6) 



By direct computations we find that there are only five combinations that transform 
homogeneously under Weyl tranformations. They are 



Of = fl<H - (d - 2) (v^V* hxT - lnT V v> Inr) , 



of 
of 



- (d-2)u a R a ^'-up 



<7 N A <7 ' 



Of 



(3.7) 

(3.8) 
(3.9) 



In the linearized hydrodynamics in flat space only the term Of contributes. For conve- 
nience and to facilitate the comparision with the Israel-Stewart theory we shall use instead 
of (|3.7| ) the term 



+ -7— r ^(V-« 
a — 1 



(3.10) 



(l/2)0 3 



20f . It is 



which, with ( |3.5| ), reduces to the linear combination: Of — Of 
straightforward to check directly that ( 3,1C| ) transforms homogeneously under Weyl transfor- 
mations. 

Thus, our final expression for the dissipative part of the stress-energy tensor, up to second 
order in derivatives, is 



if" 



-rjcr' 



+ vm 



d- 1 



+ K 



U/3 



(3.11) 



The five new constants are tr, k, Ai 5 2,3- Note that using lowest order relations n**" = —ija^, 
Eqs.(|3.5|) and Dr\ = —r]\7-u, Eq. (3.11) may be rewritten in the form 



-w^ - m 



+ 



d 



n^(v-u) 



d-l 

+ K R {fMu) - (d - 2)u a R a ^u p 



(3.12) 



This equation is, in form, similar to an equation of the Israel-Stewart theory (see Section ^). 
In the linear regime it actually coincides with the Israel-Stewart theory (6.1). We emphasize, 
however, that one cannot claim that Eq. ( 3.12; ) captures all orders in the momentum expansion 
(see Section ^). 
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Further remarks are in order. First, the k term vanishes in flat space. If one is interested 
in solving the hydrodynamic equation in flat space, then k is not needed. Nevertheless, k 
contributes to the two-point Green's function of the stress-energy tensor. We emphasize that 
the term proportional to k is not a contact term, since it contains u^. The \i,2,3 terms are 
nonlinear in velocity, so are not needed if one is looking at small perturbations (like sound 
waves). For irrotational flows A2 ; 3 are not needed. The parameter rn has dimension of time 
and can be thought of as the relaxation time. This interpretation of rn can be most clearly 
seen from Eq. ( 3. 12| ) . For further discussion, see Section ^. 



3.2 Kubo's formulas 

To relate the new kinetic coefficients with thermal correlators, first let us consider the response 
of the fluid to small and smooth metric perturbations. We shall moreover restrict ourselves to 
a particular type of perturbations which is simplest to treat using AdS/CFT correspondence. 
Namely, for dimensions d > 4 we take h xy = h xy (t, z). For d = 3, there are only two spatial 
coordinates, so we take h xy = h xy (t). Since it is a tensor perturbation the fluid remains at 
rest: T = const, = (1,0). Inserting this into Eq. ( |3.11|) we find, for d > 4, 



T xy = -Ph xy - r]h xy + rimhxy ~ \ l(d - 3)h xy + d 2 z h xy ] . (3.13) 



The linear response theory implies that the retarded Green's function in the tensor channel 
is 

G x J?' xy (uj, k) = P — irjuj + t]t u uj 2 - - \(d - 3)w 2 + k 2 } . (3.14) 
" 2 

For d = 3 there is no momentum k, and the formula becomes 

G x l xy {uj) = P - irpj + vmoo 2 , d = 3. (3.15) 

Thus the two kinetic coefficients rn and k can be found from the coefficients of the u> 2 
and k 2 terms in the low-momentum expansion of G x ^ xy {uj, k) in the case of d > 4, and just 
from the u> 2 term in the case of d = 3. 

3.3 Sound Pole 

We now turn to another way to determine rn, which is based on the position of the sound 
pole. The fact that we have two independent methods to determine rn allows us to check the 
self-consistency of the calculations. 

To obtain the dispersion relation, we consider a (conformal) hydrodynamic system in 
stationary equilibrium, that is, with fluid velocity = (1,0), homogeneous energy density 
£ = const • T d and IT"' = 0. The speed of sound is defined by c 2 = dP{e)/de. In conformal 
theory it is a constant: c 2 = l/(d — 1). Now let us slightly perturb the system and denote 
the departure from equilibrium energy density, velocity, and stress as 5s, u 1 , and IT- 7 . 



For small perturbations, one can neglect the nonlinear terms in Eq. ( 3.12 ) and the hy- 



drodynamic equations are identical to those of the Israel-Stewart theory. For completeness, 
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we rederive here the sound dispersion in this theory. To linear approximation in the pertur- 
bations, we have 



8T 00 = 5e, 5T 0i = (e + P)u\ ST^ = c 2 s 5e 5 ij + . (3.16) 

For sound waves travelling in x direction we take u x and TL XX as the only nonzero com- 
ponents of u l and IP 3 , and dependent only on x and t. Energy-momentum conservation 
implies 

d (5e) + (e + P)d x u x = , (3.17) 
(e + P)d u x + c 2 s d x (Se) + d x U xx = . (3.18) 



Eq. ( |3 .12 ) has the form 



mdoIF* + H xx = - 2 ^_*\ d x u x . (3.19) 



For a plane wave, equations ( 3.17 ), ( 3.1Sj ) and ( |3.19|) give the dispersion relation 

-uj 3 T n - iuj 2 + ujk 2 c 2 s T n + uk 2 ^ - 2 ^ + ik 2 c 2 = 0. (3.20) 

a — 1 e + P 

At small k, the two solutions of this equation corresponding to the sound wave are 

wi >2 = ±c s k - iTk 2 ± ^ (c 2 t u - ^ J A; 3 + 0(A; 4 ) , (3.21) 



where 



(3.22) 



d-l e+P 
The third solution is given by 

uj 3 = -ir~ l + 0(k 2 ) . (3.23) 

Since uj% does not vanish as k — > 0, but remains on the order of a macroscopic scale, this third 
solution lies beyond the regime of validity of hydrodynamics (see also discussion in Section ||) . 

3.4 Shear pole 

In hydrodynamics, there exists an overdamped mode describing fluid flow in a direction 
perpendicular to the velocity gradient, e.g., with u y ~ e - tujt + tkx . First-order hydrodynamics 
gives the leading-order dispersion relation, u = —ir]k 2 /(e + P). The next correction to this 
dispersion relation is proportional to k 4 and thus is beyond the reach of the second-order 
theory. This correction can be fully determined only in third-order hydrodynamics. To 
illustrate that, we shall compute this correction here, taking the second-order theory literally 
and pretending the third-order terms are not contributing. We shall than find the expected 
mismatch between this (incorrect) result and the AdS/CFT computation in the strongly 
coupled M = 4 SYM theory. 
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The perturbation corresponding to the fluid flowing in the y direction with velocity 
gradient along the x direction (shear flow) involves the variables 



u y (t,x), U xy {t,x), (3.24) 

such that we get from d^5T^ u = 

(e + P)d u y + dji xy = 0. (3.25) 



From Eq. ( 3.12 ) we find 

md U xy + U xy = -r]d x u y . (3.26) 
The dispersion relation is determined by 

-uj 2 T U -iu; + k 2 -^-^ = (3.27) 

so the shear mode dispersion relation in the limit k — > becomes 

u = -ihk 2 -ih 2 T U k 4 + 0(k 6 ), h = j^Tp- ( 3 - 28 ) 

The second solution, u = —it^ 1 + 0(k 2 ), is obviously beyond the regime of validity of the 
hydrodynamic equation (see also Section ||). 

It is easy to see that expression (|3.28| ) unjustifiably exceeds the precision of the second- 
order theory: the kept correction is 0(k 2 ) relative to the leading-order term, instead of 
being 0(k). We can trace this to Eq. ( |3.27| ), in which we keep terms to second order in u 
and k. For shear modes, however, to ~ k 2 , and the term uj 2 that we keep in Eq. ( |3.27| ) is 
of the same order of magnitude as terms 0(k A ) omitted in Eq. ( 3.27| ). The latter term can 



appear if the equation ( |3.26 ) for U xy contains a term d x u y that may appear in third-order 



hydrodynamics. This is beyond the scope of this paper. 
3.5 Bjorken Flow 

So far, we have studied only quantities involved in the linear response of the fluid, for which 
linearized hydrodynamics suffices. In order to determine the coefficients Ai2,3, one must 
consider nonlinear solutions to the hydrodynamic equations. One such solution is the Bjorken 



boost-invariant flow [16], relevant to relativistic heavy- ion collisions. 

Since hydrodynamic equations are boost-invariant, a solution with boost-invariant initial 
conditions will remain boost invariant. The motion in the Bjorken flow is a one-dimensional 
expansion, along an axis which we choose to be z, with local velocity equal to z/t. The most 
convenient are the comoving coordinates: proper time for each local element r = Vt 2 — z 2 
and rapidity £ = arctanh (z/t). In these coordinates each element is at rest: (u T ,u^,u- L ) = 
(1,0,0). 

The motion is irrotational, and thus we can only determine the coefficient Ai, but not A2 
or A 3 . 
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Since velocity is constant in the coordinates we chose, the only nontrivial equation is 
the equation for the energy density: 

De + {e + P)V • u + ^ = 0. (3.29) 

Boost invariance means that e(r) is a function of r only. The metric is given by ds 2 = 
—dr 2 + r 2 d£ 2 + dx\ and it is easy to see that the only nonzero component of V^U;, is 
V^iig = t. Using P = s/(d — 1) we can write: 

d e 



■n«. 



(3.30) 



d-1 r 

For large r, the viscous contribution on the r.h.s. in ( [3.30 ) becomes negligible and the 
asymptotics of the solution is thus given by 

d-2 



e(r) = C t + (viscous corrections), where v 



d-r 



(3.31) 



and C is the integration constant. As we shall see, the expansion parameter in ( |3,3lD is r 
Calculating the r.h.s. of Eq. ( 3.30| ) using Eq. ( |3.11| ) we find 



-rn« 



2VT]T 2 + 1v 2 I TJTTL 



2Ai 



d-3 
d^2 



+ ... . 



(3.32) 



Integrating equation ( |3.30| ), one should take into account the fact that kinetic coefficients 
rj, Tn and Ai in Eq. ( |3.32| ) are functions of s, which in a conformal theory are given by the 
following power laws: 



v = Cm ( ^ ) 



e\(d-i)/d 



C 



-l/d 



Q / £ \(d-2)/d 



Ai = CX", ( - ) 



(3.33) 



where, for convenience, we defined constants 770, and A?, and we chose the constant C to 
be the same as in Eq. ( |3.31| ). Integrating Eq. ( p.30 ) we thus find 



£(T) 

C 



,-2+v 



2% t~ 2 + 



2(d-l) 2 



d-2 



2A 



d-3 
d^2 



r- 2 -^ + . . . . 



(3.34) 



d ,u d-1 

In Section 4.4 we shall match the Bjorken flow solution in the strongly-coupled N = 4 
SYM theory found in [|l7j (see also [p^]) using AdS/CFT correspondence and determine Ai 
in this theory. 

In order to compare our results to the ones obtained in Ref. [17], we shall write here the 
equations of second-order hydrodynamics using also the alternative representation ( 3.1 2| ) for 
11^ in ( 3.30 ). We obtain the following system of equations for the energy density and the 

-H| (see fn§; c.f. @ for Ai = 0): 



component of the viscous flow, which we define as $ 



d e + £ 

d— 1 T T 



d— 3 Ai l9 
i <|) z 

d-2 i] 2 



(3.35) 
(3.36) 



d-1 r d- 1 t 

As should be expected, the asymptotics of the solution of this system coincides with Eq. ( p.34| ). 
Equation ( 3.36 ) is different from the one used in 17] by the last two terms proportional to 
rn and Ai. 
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4. Second-order hydrodynamics for strongly coupled M = 4 supersymmetric 
Yang-Mills plasma 



In this Section, we compute the parameters rn, n, and Ai of the second-order hydrodynamics 
for a theory whose gravity dual is well-known: J\f = 4 SU(iV c ) supersymmetric Yang-Mills 



theory in the limit N c — > oo, g 2 N c — » oo |2l], 22 1. According to the gauge/gravity duality 
conjecture, in this limit the theory at finite temperature T has an effective description in 
terms of the AdS-Schwarzschild gravitational background with metric 

ir 2 T 2 L 2 L 2 
dsl = {-f(u)dt 2 + dx 2 + dy 2 + dz 2 ) + - - - du 2 , (4.1) 

U 4:J{U)U Z 



where f(u) = 1 — u , and L is the AdS curvature scale [14j. The duality allows one to compute 
the retarded correlation functions of the gauge-invariant operators at finite temperature. The 
result of such a computation would in principle be exact in the full microscopic theory (in the 
limit N c — > oo, g 2 N c — > oo). As we are interested in the hydrodynamic limit of the theory, 
here we compute the correlators in the form of low-frequency, long-wavelength expansions. 
In momentum space, the dimensionless expansion parameters are 

"=2Sr <<1 ' « = ir <<h < 42) 

Comparing these expansions to the predictions of the second-order hydrodynamics obtained 
in Sections |3.2|, 3.2 and 3.5 for d = 4, we can read off the coefficients rn, k, Ai. 



One must be aware that the M = 4 SU(iV c ) supersymmetric Yang-Mills theory posesses 
conserved ii-charges, corresponding to SO(6) global symmetry. Therefore, complete hydro- 
dynamics of this theory must involve additional hydrodynamic degrees of freedom - i?-charge 
densities. Our discussion of generic conformal hydrodynamics without conserved charges can 
be, of course, generalized to this case. This is beyond the scope of this paper. Here we 
only need to observe that since the i2-charge densities are not singlets under the SO (6) they 
cannot contribute at linear order to the equations for T^ v . These contributions are therefore 
irrelevant for the linearized hydrodynamics we consider in Sections [T^, |3.3| and |3.4j . For the 
discussion of the Bjorken flow in Section |3.5| they are also irrelevant, since (and as long as) 
we consider solutions with zero i?-charge density. 

4.1 Scalar channel 

We start by computing the low-momentum expansion of the correlator G^ y (u,k). To 
leading order in momentum, this correlation function has been previously computed from 



gravity in [23, 24]. Following [p4]j , here we obtain the next to leading order term in the 
expansion. 

The relevant fluctuation of the background metric ( [4.1| ) is the component (f> = h% of the 
graviton. The retarded correlator in momentum space is determined by the on-shell boundary 
action 

S tot [H , k) = run(sZ V ndary [Ho, e, k] + S c , t .[H , e,k]\, (4.3) 
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following the prescription formulated in [£3|. Here Ho(k) = H(e,k) is the boundary value 
(more precisely, the value at the cutoff u = e — > 0) of the solution to the graviton's equation 



of motion (Eq. (6.6) in \2§) 



H(u,k) = H (k) 



<Pk(u) 



(4.4) 



A perturbative solution (j)k{u) to order j» 2 , q 2 is given by Eq. (6.8) in 24]. The gravitational 
action (Eq. (6.4) in j24|) reduces to the sum of two terms, the horizon contribution and the 
boundary contribution. The horizon contribution should be discarded, as explained in 23] 
and later justified in 25]. The remaining boundary term, S^^^ylHo, e, k], is divergent 
in the limit e — » 0, and should be supplemented by the counterterm action S c .t.[Ho,e,k] 
following a procedure known as the holographic renormalization. 1 In the case of gravitational 
fluctuations, the counterterm action is [27] 



S, 



ct 



3iV 2 
"4vr 2 L 4 



dV=7( l + y p ~ S ( pklPkl loge ' ' 



(4.5) 



where jij is the metric ( |4. ip restricted to u = e, and 



Pi, 



P = l ij Pii , 




Evaluating (4.5), we find the total boundary action 2 



(4.6) 



5, 



tot 



7T 2 iV 2 T 4 



Va 



H{e)H'{e) H 2 {e) (q 2 - xo 2 )H 2 (e) 



+ 



+ O(t» 3 ,toq 2 ) + O(e). (4.7) 



The boundary action ( ]4.7| ) is finite in the limit e — > 0. Its fluctuation-independent part is 

en 



o 

tot 



-PVa, where P = 7r 2 iV 2 T 4 /8 is the pressure in J\f = 4 SYM, V4 is the four-volume. 
The part quadratic in fluctuations gives the two-point function. Substituting the solution 



Q4.4D into Eq. Q4.7]) and using the recipe of IpSj] , we find 

„2 AT2rri 



7T JVgT 4 



ixx> - t» 2 + q 2 + to 2 In 2 



+ 0(to 3 ,t»q 2 ) . 



(4.8) 



Comparing Eq. ( f4.8j ) to the hydrodynamic result ( |3.14| ) we obtain the pressure |28|] , the 
viscosity p9| and the two parameters of the second-order hydrodynamics for J\f = 4 SYM: 



P= ^T\ v = lN 2 T 3 , r n 



In 2 



2vrT 



ttT 



(4.9) 



lr rhe holographic renormalization corresponds to the usual renormalization of the composite operators 
in the dual CFT. 

2 Terms quadratic in H in Eq. ( |4.7| ) should be understood as products H(—u>, —k)H(u>, k), and an integration 
over ui and q is implied. 
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4.2 Shear channel 



The dispersion relation ( 3.2S| ) manifests itself as a pole in the retarded Green's functions Gj^ ty , 
G$j XV , Gxy,xy i n the hydrodynamic approximation. To quadratic order in k this dispersion 



relation was computed from dual gravity in Section 6.2 of Ref. [24]. Here we extend that 
calculation to quartic order in k. This amounts to solving the differential equation for the 
gravitational fluctuation G(u) [24] 



G" 



(2u 

V7 



G' + -l2 + 



nv 



q 2 t» 2 [4-u(l+u) 2 ] 



+ 



4uf 



G = 



(4.10) 



q 2 . The solution G(u) is supposed to be regular at 



l-uj- ■ f 

perturbatively in t» and q assuming to 
u = 1 p4|] . Such a solution is readily found by writing 

G(u) = G (u) + wGi(u) + q 2 G 2 (u) + to 2 G 3 (u) + n>q 2 G 4 (u) + q 4 G 5 (u) + • • • (4.11) 

and computing the functions Gi(u) perturbatively 3 . The functions Gi(u) are given explicitly 
in Appendix [A]. To obtain the dispersion relation, one has to substitute the solution G(u) 
into the equation (6.13b) of [24] and take the limit u — > 0. The resulting equation for to, 

Aim 



q 4 + 2q 2 



iwq 2 In 2 + 2» 2 ln2 



0. 



(4.12) 



has two solutions one of which is incompatible with the assumption to *C 1. The second 
solution is 

i(l -ln2)q 4 



!5I 

2 



+ 0(q b 



If we naively compare Eqs. (|3.28| ), ( 4.13 ), we would get rn = (1 
inconsistent with the value obtained from the Kubo's formula, Eq. (14. 



(4.13) 

ln2)/(2vrT), which is 



As explained in 

Section 3.4, this happens because the 0(k i ) term in the shear dispersion relation is fully 
captured only in third-order hydrodynamics. In other words, we confirm that Eq. ( 3.28Q has 
an error at order 0(k A ). 

4.3 Sound channel 

The sound wave dispersion relations ( p.21 ) appear as poles in the correlators of the diagonal 
components of the stress-energy tensor in the hydrodynamic approximation. These correla- 
tors and the dispersion relation to quadratic order in spatial momentum were first computed 
from gravity in |30[| . A convenient method of studying the sound channel correlators was 
introduced in In this approach, the hydrodynamic dispersion relation emerges as the 

lowest quasinormal frequency of a gauge-invariant gravitational perturbation of the back- 



ground (4.1). According to the sound wave pole is determined by solving the differential 
equation 



Z" - 



+ 



3t» 2 (l 



u 



+ q 2 (2u 2 -3-u 4 -3) 



n/(3to 2 + q 2 (u 2 -3)) 
3ro 4 + q 4 (3 - Au 2 + n 4 ) + q 2 (4n 5 



Z' 



4u 3 + 4u 2 i» 2 - 6t» 2 



n/ 2 (3tD 2 +q 2 (n 2 -3)) 



Z = 



(4.14) 



3 Note that, for u real, G*{u, —to) = G(u, t»). This implies ImGo,2,3,5 = 0, ReGi,4 = 0. 
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with the incoming wave boundary condition at the horizon (u = 1) and Dirichlet boundary 
condition Z(0) = at the boundary u = 0, and taking the lowest frequency in the resulting 



quasinormal spectrum. The exponents of the equation ( 4.14 ) at u = 1 are ±irx>/2. The 
incoming wave boundary condition is implemented by choosing the exponent — ivo/2 and 
writing 

Z(u) = r m/2 X{u) , (4.15) 
where X(u) is regular at u = 1. Thus we obtain the following differential equation for X(u) 



X" + 



+ 



2uixo 1 + u 
(1 + u + u 2 )vo 2 



4q 2 u 



3u> 2 + q 2 (u 2 - 3) 

,2„3 



X' 



4q 2 u 6 {l + ixo) 



u{l + u)f uf u/(3» 2 +q 2 (u 2 -3)) 



\X = 0. 



(4.16) 



This equation can be solved perturbatively in to <C 1, C) <C 1 assuming vo ~ q (the expected 
scaling in the sound wave dispersion relation). Rescaling vo — > Ato, q — > Aq, where A C 1, we 
look for a solution in the form 



X(n) = X (u) + AXi(ii) + A 2 X 2 (n) + • • • . 



(4.17) 



The functions Xj(u) are written explicitly in Appendix |A[ The Dirichlet condition ^(0) = 
leads to the equation for ro(q): 



• 2 , <T 
't»q + y 

ro 2 q 2 , 2 



^ + ^( vr 2 _121n 2 2 + 241n2) -^-(21n2 
2 16 v ; 12 v 



48 



vr 2 - 12 In 2 2 + 48 In 2) = 0. 



(4.18) 



To order q 3 , the solution is given by 



(4.19) 



This is the dispersion relation for the sound waves to order q 3 . The complete dispersion 
relation can be obtained by solving the equation ( [4.14 ) numerically pl| . The sound dispersion 
curve is shown in Fig. [p. Comparing Eq. ( [4.19| ) to Eq. ( |3.21 ) we find the relaxation time Tn 
for the strongly coupled J\f = 4 SYM plasma: 



2 - In 2 
2ttT 



(4.20) 



The result ( 4.20)) coincides with the one obtained in Section 4.1 , which is a nontrivial check 
of our approach. 
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Figure 1: Sound dispersion c s — c s (q) in Af = 4 SYM plasma. The dark (blue) curve shows the sound 
speed dependence on wavevector, c s (q) = Ret»/q, with c s (0) = l/VS (this plot is based on numerical 
data first obtained in |3lJ] ) . The light (red) curve corresponds to analytic approximation derived from 



Eq. (4.19) and valid for sufficiently small q. 



4.4 Bjorken flow 

In order to determine Ai, we match Eq. (3.34) with the solution found by Heller and Janik (l 
given by 4 



e(r) 



N 2 

± "c 
2rr 2 



^ /3 - 2 "^ 2 + ^ /3 'T"»+ 36^3 



10 



6m2- 17 



with 7/n 



1 



^/2 3 3 / 4 



(4.21) 



Matching by using C = N 2 /(2ir 2 ), and r n = (2 - ln2)/(27rT) from Eq. ( [4.20D , together with 
e = 3vr 2 iV c 2 r 4 /8 and Eq. <^3^) gives 

x _ V 



2vrT 



(4.22) 



Note that Heller and Janik jl7|] found a different value for rn since they matched to the Israel- 
Stewart equations for hydrodynamics, and not the more general (nonlinear) equation ( |3.12 ). 



5. Kinetic theory 

Our analysis should be valid not only for the strongly coupled M = 4 SYM theory, but also 
for all theories with conformal symmetry. In particular, it should be valid also for weakly 
coupled CFT like the SYM theory at small 't Hooft coupling, or QCD at sufficiently large 
Nf at the Banks-Zaks fixed point jj^]. In these cases, one expects that it is possible to 
understand and compute the second-order transport coefficient from kinetic theory. We set 
d = 4 in this Section. 



4 The quantities in Eq. (4.21) can be thought of as dimensionless combinations of quantities in Eq. (3.34) 
with an appropriate power of an arbitrary scale parameter to: t/to, stq, tiqtq , Ctq v etc. Due to conformal 
invariance, a rescaled solution is also a solution, and the scale to can be used instead of the integration constant 
C, to parameterize the solutions in Eq. (3.34). 
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5.1 Setup 



Since we are to discuss conformal transformations, our starting point is the classical Boltz- 
mann equation in curved rather than flat space-time 



P dx n ^P P dp \ 



f(p,x) = -C[f], 



(5.1) 



where f(p,x) is the one-particle distribution function, p^ is the particle momentum, Y^ u are 
the Christoffel symbols and C is the collision integral. One can easily show that conformal 
transformations are a symmetry of the Boltzmann equation if particles are massless (p 2 = 
P^Pti = 0) and the collision integral transforms as C[f] — > e 2ui ( x >C[f]. 

Hydrodynamic equations are obtained by taking moments with respect to the particle 
momentum p^ of Eq. ( |5.1| ). More precisely, acting with J d\ = j d 4 p5(—p 2 )9(p°), where 9 is 
the step-function, on Eq. (fOI) one obtains 



dxV^g 



d 
dp x 



dx^gC[f] 



which upon partial integration leads 34] to 



V, 



dxV=gC[f}. 



(5.2) 



(5.3) 



We recall here that is the (geometric) covariant derivative. In theories with conserved 
charges or if only elastic collisions are considered, f dxC[f] = and Eq. ( p^ ) becomes the 
conservation of the particle current in theories with conserved charges. Conservation of the 
energy-momentum tensor 5 

T^= [ dxpVV=0/(p,aO (5-4) 



follows from Eq. ( |5.1| ) upon action of J dxp u and the requirement J dx\/ —gp u C [/] = 0, 



(5.5) 



Acting with J dxp v p k on Eq. ( |5.1[ ) gives the first equation with non-trivial contribution from 

(5.6) 



the collision integral [36], 



VfiX ^ x = I 



where 



dxp^p u p x \^gf{p,x) , 



dxp v p X y 



-gc[f}. 



(5.7) 
(5.8) 



5 Note that sometimes p M is traded by the introduction of a "local momentum" j35| and as a consequence 
T^ v would be defined without a factor of yj— g and the form of the Boltzmann equation (5.1) changes. 
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Similarly, an infinity of higher moment equations of the form 



V fl X' wlvava " = I uiv * va ~ (5.9) 

also follow from Eq. (|5.1|) . 

Splitting the out-of-equilibrium particle distribution function into an equilibrium and 
non-equilibrium part 

fiP, X) = /eq(p, X) (1 + 5f(p, x)) , (5.10) 

one defines an equilibrium energy-momentum tensor 

T% =T^ = J dxp»p u V=gf cq (p, x) , (5.11) 

and a non-equilibrium component IP 41 ' = T^ v — Teq , which by construction is both sym- 
metric and traceless. We shall assume that the equilibrium distribution function f eq (p,x) = 
feq(—u(x) -p/T(x)) depends on local temperature and velocity T, u M , which are defined such 
that the equilibrium distribution has the same energy and momentum density as / in the rest 
frame defined by u^, 

J d X V=9P»{u vP v ) (f - /eq) = 0. (5.12) 

This implies that u^W v = 0. 
5.2 Moment approximation 

While the full hierachy of moment equations should correspond to the original Boltzmann 
equation, it is too complicated to be treated exactly. However, an approximate evolution 
equation for systems not too far from equilibrium may be constructed. The approximation is 
similar to the Grad's 14-moment method |3^| . 
We decompose Sf into spherical harmonics, 

oo 

^/ = E4?... w (0p mi ---p m s e = -^, (5.i3) 

1=0 

where f^...^^) are fully symmetric, orthogonal to n M , and traceless over any pair of indices. 
By construction, the I = 0, 1 parts satisfy the constraints Eq. ( 5.12j ). The approximation is 
now specified by the following assumptions (c.f. [j38f): 

• the system is sufficiently close to equilibrium that the collision term is linear in 5f 

• all contributions / > 2 are sub dominant 

• for / < 2 and expanding in some basis, all £ dependent terms are subdominant. 



- 19 - 



This implies that 
and 



5/( P ,x)~T-yyn^ + 0(n 2 ), 



(5.14) 



T <ux> ^ T 2( x ) n <^A> (x) + {Tl 2 ), (5.15) 

where subdominant terms have been labelled as 0(II 2 ). It would be interesting to use nu- 
merical techniques such as in Ref.|H |C| to test the correctness of Eq. |TlJ). 
Splitting X^ x into an equilibrium and non-equilibrium part, one finds 



X cq X 



where perm, denotes all non-trivial permutations of indices, and 



rpt) 



u> 1 u ,J u x + const x (A^u x + 



■ perm. 



where {^1^2 ■ ■ ■ Hn) denotes symmetrization with respect to the indices m,fJ,2, 
jection <> on the moment equation ( |5.6| ) thus gives 



, (5.16) 

(5.17) 
,fj, n . Pro- 



n^ A + r n 



U uX DlnT + AXA^£>n Q/3 + Ii uX V ^ + 2W <U V 'u x> 



r/a uX + 0(U 2 ), (5.18) 



where the proportionality constants have been denoted by r/ and rn, respectively (the ratio of 
these can be calculated when specifying / eq , c.f.[^]). Introducing the completely symmetric 
tensor 



one can decompose 



m <v V^u x> 



-u a( - u n x l + ir^e x) - -iL a ^A vX e. 



a/3- 



Rewriting 
such that 

we find 
U uX = -rja 

where DlnT 



1 T\a<v-rx\> 

n» <u v, lU x> = -n a ( u n x l + -ir A v 7 u 7 + orn 3 

3 2r] 



DU <uX> + -U uX (V -u) 



+ 2r n n Q (^ A ) + ^U a<u U x> + 0(U 3 



(5.19) 

(5.20) 
(5.21) 
(5.22) 

(5.23) 



4(V • u) + C(n 2 ) has been used. 



Eq. ( 5.23j) , which was derived from kinetic theory here, corresponds to the more general 
Eq. ( |3.12[) with A2 = —2tui] and A3 = K = 0. Note that Ai contains a contribution from 
Eq. (5.22) as well as from the collision integral Eq. ( |5,15| ) (see below). What is commonly 
referred to as Israel-Stewart theory amounts to setting Ai = 0. Most of the time, also the 
terms involving V • u and the vorticity W are dropped. However, note that simply dropping 
terms involving V • u ruins the conformal symmetry of the equation, and thus the resulting 
equation cannot be the correct hydrodynamic description of the system dynamics beyond 
leading order. 
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5.3 The structure of the collision integral 

In this subsection we study the structure of the collision integral Eq. ( |5.15 ) for a simplified 



£ £ 

model where C = (u ■ p) - — — . We will use a gradient expansion similar to the Chapman- 



Enskog method (c.f. Q). 
Let us decompose / into 

/ = /eq(-« • P/T) (1 + /l + h + ■ ■ ■) , (5-24) 

where /i,/2 represent terms of first and second order in gradients, respectively. Solving 
Eq. iteratively in gradients we find 

fi = — -j^P^P V M — , 

P-Ufeq T 



r 2 f ?/ 2r 2 f 



p^p a V u V^ + ^ipf^pav (^) V^lnT. (5.25) 



(p- U ) 2 / e / T (p-n) 2 / eq " r ' "VT 

From Eq. ( |5 . 1 5| ) and conformal symmetry, to second order in gradients the collision integral 
j<-/S> can con tain terms <7^ 7 <7 5>A and Da <lX> + |cr' 1 '^(V • u) but (in particular) not J7 7<5 or 
i? 7<5 since these terms would involve anti-symmetrization of indices which is not allowed by 
Eq. flpp . 



This indicates that the terms involving k, A3 in Eq. ( 3.12 ) are not contained in the Boltz- 
mann equation. The Boltzmann equation is only an approximation of the underlying quantum 
field theory, so it is possible that these terms - which are second order in gradients - have 
been lost in this coarse-graining process. It may be possible to compute the coefficients of 
these terms for QCD in the weak-coupling regime by going beyond the lowest order gradient 
expansion given in p^ ]. 

6. Analysis of the Miiller-Israel-Stewart theory 
6.1 Causality in first order hydrodynamics 

It is instructive to compare the second-order conformal hydrodynamics to the Miiller-Israel- 
Stewart theory. Miiller 0] and independently later Israel and Stewart H ^, 1C], considered 



how to extend the 1st order hydrodynamics. Their primary motivation was to eliminate 
the apparent relativistic acausality of the 1st order hydrodynamic equations. Formally, the 
acausality is the result of the fact that the 1st order hydrodynamic equations are not hyper- 
bolic [43, 10, The problem is most clearly seen by considering the linearized equation for 



a diffusive mode (e.g., shear stress or charge diffusion), which is first order in temporal but 
second in spatial derivatives. A discontinuity in initial conditions for such a mode propagates 
at infinite speed. In other words, the influence of an initial condition at a point in space is 
instanteneously felt by any other point. 
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It should be clear, however, as emphasized, e.g., by Geroch |^5|, [|(| and others JI?]] that 
the modes which defy causality are those which are not supposed to be described by hy- 
drodynamics (i.e., microscopically short wavelengths, which is clear when one thinks about 
discontinuities). Nevertheless, for numerical simulations of relativistic hydrodynamic systems 
such superluminal propagation is a nuisance because in such simulations one extrapolates 
hydrodynamic equations to the microscopic scale, even though the modes, or the configura- 
tions, which are being studied are hydrodynamic. For example, superluminal propagation 
makes posing initial value problem difficult: even if the initial hypersurface is space-like, the 
initial values at different points can influence each other and an attempt to specify them 
independently leads to unacceptable singular solutions |4^, f47| . 

Since the problem lies in the domain where the theory is not applicable, one can safely 
modify the theory in this domain, without disturbing physical predictions. This is the essence 
of the solution which Miiller and Israel proposed by extending the set of variables. The 
resulting system of equations is hyperbolic. Here we shall write down explicitly the system 
of equations of Israel and Stewart, restricting to the case of conformally invariant system 
without a conserved charge that we study in this paper. 



6.2 Hydrodynamic variables and second order hydrodynamics 



As we have already emphasized in Section |2.2| the hydrodynamics should be viewed as a 
controllable expansion in gradients of the hydrodynamic variables. The choice of the variables, 
or fields, can be aided by applying the requirement that a linearized system of equations has 
solutions whose frequency vanishes in the hydrodynamic limit, i.e., when the wave vector k 
vanishes. We call such linearized modes the hydrodynamic modes. Fluctuations of conserved 
densities are automatically hydrodynamic because their equations are conservation laws and 
constant fields (u = 0, k = 0) are trivial solutions of them. 

Hence, for a system without conserved charges the set of hydrodynamic variables consists 
of the densitites of energy and momentum, represented by 4 independent covariant variables 
e and (u-u = —1). All other quantities in hydrodynamic description are instantaneous 
functions of these variables and their derivatives, such as, e.g., IL^ (Section p.2| ). 

How should one extend 1st order hydrodynamics to higher derivatives? The systematic 
way, as we argued in Section 2.2 and|3|, is to continue the expansion (2.16) and add all possible 



terms of the second order in derivatives, as we did in Eq. (3.11). 

Instead, Miiller, Israel and Stewart take a more phenomenological point of view. They 
consider Ii^ u — the viscous part of the the momentum flow - as a set of independent additional 
variables. The equations for these variables are not given by any exact conservation laws, 
but by phenomenological expansions in the set of independent variables, which now includes 
also m v : 

tjlDW = -BT - rja^ . (6.1) 

The first term in Eq. ( |6.lD has a simple intuitive meaning: in the absence of velocity gradients 
(a^ v = 0) the viscous momentum flows H fJ,u do not vanish instanteneously (as in Eq. ( 2.16[ )), 
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but relax to zero on a microscopic but finite timescale rri- The 5 equations (6.1) together 
with 4 conservation laws V^T^ = form the system of Miiller- Israel-Stewart equations for 9 
variables: s, u^ 1 and Ii^ v . (For a non-conformal system with a conserved charge this number 
becomes 14.) 

In the phenomenological laws in Eq. ( |6.1| ) one usually considers only terms linear in 
the variables and u^. There is a priory no reason to neglect nonlinear terms. By 
comparing Eq. fl6.1| ) with Eq. fl3.12| ) we see that the conformal invariance requires presence 
of terms proportional to U^ u (V-u), which are non-linear, but contain the same number of 
derivatives. These terms are beyond the standard linear Israel-Stewart phenomenological 
theory. In addition, bilinear terms proportional to Aj are also allowed to the same order 
in derivatives. Such terms are relevant for simulations of the strongly coupled quark- gluon 
plasma in heavy ion collisions. 

The term proportional to k, which vanishes in flat space, has not been considered by 
Israel and Stewart but, as we have seen, is necessary to determine the correlation functions 
of stress-energy tensor. 

Note that in this scheme both IP y and a^ u are of the same, i.e., first order in the 
expansion around equilibrium. The term DW LU contains one more derivative compared to 
U^ u and is thus of the second order. Without loss of precision, to second order, one can trade 
DU^ U for —D(r]a tJ " v ) or vice versa. Similar substitutions can be made in other second-order 



terms we found, as we did when going from Eq. ( 3.11| ) to Eq. ( |3.12|) . Therefore, within their 



precision, equations of Israel-Stewart ([O]) (or, in general nonlinear case, Eq. ( |3.12D ) give the 
same result as the systematic expansion in derivatives. 

6.3 Causality and the domain of applicability 

The attractive feature of introducing new variables is that the resulting equations are now 
first order in derivatives and, most importantly, they are hyperbolic. This means that dis- 
continuities propagate with finite velocities even in the shear channel. For the shear channel 
this velocity (i.e., the characteristic velocity fl43| , 49, Q) can be easily obtained from the 
dispersion relation ( 3.27| ) by taking k — > oo: 



Vdisc = ^I^Tp-y (6 ' 2) 



Although the Israel-Stewart system of equations fl6.1|) or our equations ( 3.12 ), have at- 
tractive features from the point of view of the mathematical formulation, and are especially 
suitable to, e.g., numerical simulations, care should be taken attributing physical significance 
to this fact. The domain of applicability of these equations is still the hydrodynamic do- 
main: u, k must be small compared to microscopic scales. The second order hydrodynamic 
equations increase the precision compared with the first order equations, but only if we stay 
within the hydrodynamic domain. 

In practice, it is convenient to use equations which are mathematically well-behaved 
even where they lose physical significance. However, care should be taken when examining 
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the solutions by always considering only their features in hydro dynamic domain - slow and 
long-wavelength modes. In particular, the velocity in Eq. ( |6.2[) does not correspond to any 
physical propagation. Similarly, the super luminal propagation which one recovers according 
to Eq. ( |6.2[ ) in the first order theory when rn — ► is the result of extrapolating the theory 
outside the hydrodynamic domain. 

Nevertheless it is worthwhile to note that, with the value of rn in strongly coupled J\f = 4 



SYM that we find in Eq. (4.9), the characteristic velocity ( |6.2j ) equals l/y / 2(2 — log 2) 



0.6. . ., i.e., less than the velocity of light. Therefore, the system of second order equations 
we wrote down can be used in, e.g., numerical simulations without additional modifications 
often needed to ensure relativistic causality and prevent occurence of singular solutions. 

6.4 Entropy and the second law of thermodynamics 

Let us consider the question of how the second law of thermodynamics is obeyed by the 
second order hydrodynamics. For that purpose take the projection of the energy-momentum 
conservation equation on u u : 

= -u v V^ v = De + (e + P)V-u + II^V^, (6.3) 



where we used definition Eq. (3.1), u-u = — 1 and u v Ii^ u = 0. For a system without a 
conserved charge, the thermodynamic entropy density s is a function of the energy density 
such that ds = de/T, and it also obeys sT = e + P. Thus, Eq. (|6.3|) can be writen as 



TVJsu^) = -n^W. (6.4) 



Since s is the entropy in the local rest frame, equation ( |6.4[ ) expresses, in a covariant form, 
the rate of entropy production in the local rest frame. 

For a conformal system the tensor Yi^ is traceless and one can replace 'V^u u on the r.h.s. 



of Eq. (6^4) with /2. Using the first order hydrodynamic relation ( |3.2|) one then finds 



V^(su^) = ^ a^a^ + (3rd order terms). (6.5) 



Thus, if r) > 0, the entropy increases, provided the 3rd order terms on the r.h.s. of Eq. fl6.5|) 
are negligible compared to the 2nd order term written out. This is always true within the 
domain of validity of hydrodynamics. 

Miiller and Israel observed |7], [8| that the third order terms in Eq. ( |6.5[ ) in their theory 
can be written as the divergence of a current. Indeed, even a complete, conformally covariant, 
term proportional to rn in Eq. (|3.12j ) can be written in such a way. Solving (|3.12j ) for a^ u 
and substituting into Eq. ( |6.4[) we find 



1 + A 1 Of + A 2 Of + A 3 On + ... > (6-6) 



2rjT 



- 24 - 



where we used rn/?? = const-T _rf and the lowest order relation D In T = — (d — l)V-u. The el- 
lipsis in Eq. ( |6.6| ) denotes 4-th order corrections. Therefore, defining non-equillibrium entropy 
as 

Snoneq = 8- ^f^-afi^-^ (6.7) 

one can cancel the 3rd order term proportional to rn in V At (s non eqW At )- The correction to the 
equilibrium entropy in Eq. ( |6.7[ ) has an intuititive meaning - a non-homogeneous state of the 
system, in which ^ 0, has smaller entropy than the equilibrium state. 

The remaining terms, such as e.g., kIL^O^ 1 ' / (j]T) , do not appear to be total derivatives. 
They are also not positive definite. However, this fact cannot be used to conclude that, e.g., 
k must be zero. Our explicit AdS/CFT calculation shows that As we discussed above, 

the 3rd order terms in Eq. ( |6.5D do not violate the second law of thermodynamics if we stay 
within the domain of applicability of hydrodynamics. In this domain the 3rd order terms 
must be small compared to the second order term on the r.h.s. of Eq. Q6.5| ), which is positive 
definite. 

Further detailed discussions on the issue of the local entropy current can be found in 5C, 

H- 

6.5 Additional non- hydro dynamic modes 

Another interesting consequence of introducing more variables, a la Miiller- Israel-Stewart, 
is that the number of modes, or branches of the dispersion relation uj(k) increases, as we 



have seen in Sections [3^ and 3.4. As should be expected, the additional poles are not 
hydrodynamic: those frequencies oj(k) do not vanish as k — > 0, but remain on the order of 
the microscopic scale. It should be clear from the discussion above that the position of these 
poles need not be predicted correctly by the second-order theory - they lie outside of the 
regime of its validity. 

In fact, now with the knowledge of the position of Green's function singularities in J\f = 4 
SYM at strong coupling |3^] we can say that there are infinitely many such poles. They are 



given by the solutions of equations such as ( [4.10 ) or (|4.14 ). Only the lowest branch u>(k) 



can be matched by hydrodynamic theory. To describe correctly the full Green's function one 
needs to introduce infinitely many degrees of freedom - to describe infinitely many poles. Any 
theory of finite number of degrees of freedom is a truncation. This truncation is controllable 
only for the hydrodynamic variables, which describe the poles with frequencies vanishing as 
k — > 0. The controlling parameter is the ratio of these frequencies to a microscopic scale, i.e., 
T in the conformal theory, and the precision can be, in principle, increased by increasing the 
order of the expansion in this parameter. 

Conceptually, let us imagine that we did succeed in writing the infinite set of extended 
hydrodynamic equations for infinitely many variables, mentioned in the previous paragraph. 
It is easy to realize that in a theory with gravity dual this set will be mathematically equivalent 
(in the linear regime) to differential equations ( 4.10) ) or ( [4.14| ). The set of infinitely many 



4-dimensional fields is represented by a 5-dimensional field in these equations. 
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7. Conclusion 

We have determined the most general form of relativistic viscous hydrodynamics of a confor- 
mal fluid (with no conserved charges) to second order in gradients. We find that conformal 
invariance reduces the number of allowed terms relative to more general, non-conformal, hy- 
drodynamics. As already known, at first order in gradients only one kinetic coefficient, the 
shear viscosity rj, enters the equations. At second order we find five allowed terms with 
coefficients rn (customarily referred to as relaxation time), k, X±, X2 and A3. 

The general viscous hydrodynamic equations we obtained can be matched to AdS/CFT 
calculations in strongly coupled N = 4 supersymmetric Yang-Mills theory, and for this theory 
we thus determined three of the five second-order coefficients: rn, k and Ai. We also find that 
for a weakly coupled conformal plasma describable by the Boltzmann equation, two of the 
coefficients vanish. However, at least one of these coefficients, i.e., k, is not zero for strongly 
coupled N = 4 super Yang-Mills theory. It would be interesting to understand how this 
coefficient emerges as the approximation of the Boltzmann equation breaks down at large 
coupling. 

We emphasized the already known fact that the equations of the Muller-Israel-Stewart 
theory, despite their appearance, are only applicable in the hydrodynamic regime, where their 
predictions coincide with those of the second-order gradient expansion. We also pointed out 
that variants of the Miiller- Israel-Stewart theory used in numerical simulations of relativistic 
plasmas frequently miss terms which are not only allowed, but also required for conformally 
invariant theories. If the quark-gluon plasma is approximately conformal, then the second- 
order hydrodynamic equation found in this paper should be used instead. One may hope that 
the values of the kinetic coefficients rn and Aj, found in N = 4 SYM theory, may serve as 
crude estimates for their values in the strongly coupled regime of the quark-gluon plasma. 
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A. Perturbative solutions of the shear and the sound mode equations 



The shear mode 

The functions Gi(u) entering the perturbative solution ( 4.1 1| ) of the equation ( 4.1C| ) are 



G 2 (u) 



C(l-u) 



(A.l) 



G 3 (u) = - —I 6vr 2 u - 24(u + 1) In 2 - iUiru In 2 - 6n In 2 2 + 18u In 2 (14 - 1) 
48 \ 

+ 24(i4+l)ln(u + l) + 12i4ln21n-!-^ - 12nln(l + u) In + 6u i n 2 2_±J£ 

1 — U 1 — 14 1 — 14 



- 24uLi 2 



1-14 



+ 12i4 ln(n - 1) In 2 - 2 ln(l - u) - m 



GAu) = — ( -4vrii-4i(l + 3i4)ln2 + 4iln(l + i4) + 16iuln^-^ + 2mln(n-l) (in — — 
16 \ u \ 1 — u 



1+14 1+14 

Aiu In h 2iru In ■ 



1 - 14 



G 5 (u) = jll-u-2uln 



1-14 

l + u\ 



. 1 + U , n 1 + U 

2iu ln(l + it) In h 2ii4 In 2 

1 — u 1 — 14 



4iiiln(i4 — 1) 



(A.2) 



where C is a constant, Li2(z) is a polylogarithm. 

An alternative way to obtain the dispersion relation ( |4.12| ) is the following: the functions 
Gi(u), i = 0, 1, ..5 satisfy the inhomogeneous differential equations 



1 - u^Gi - 2uG i + 2d = Fi(u) 



(A.3) 



with Fq = 0, F% = — i(l + u)G' — i/2Go, etc. The homogeneous part of ( |A,3j ) is the Legendre 
differental equation with the Legendre functions Pi{u) = u and Qi(u) = § In — 1 as 



solutions. Therefore Go = Ci4, and for i > 1 



Gi{u) = Px(«) f Q^Fi^du' - Qi(u) / Pi(u')Fi(u')du' , 

J u J u 

regular at u = 1. Finally, the values at u = are obtained by 



(A.4) 



Gi(u = 0) = / uFi{u)du , 
Jo 



i.e. G (0) = 0,Gi(0) = -iC, and 



(A.5) 



G 2 (0) = / G (u)du = C/2, 
J o 



GM^f 



(2 + 3u) In i^H + 7n — 

2 l + ii 



dn2 



tin = C , (A.6) 
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etc., and hence we find Eq. ( [Hp . 
The sound mode 

The functions Xi(u) of the perturbative solution ( 4.17 ) of the equation ([4.161 ) are 



X (u) 



(q 2 +q 2 « 2 -3t» 2 )C 
4q 2 



iCwf(u) 



(A.7) 



Xo(u) 



C 



2q 2 8 - 8u - mt(1 + u 2 ) - (1 + u 2 ) 2 In 2 



48q 2 

+ 3t» 4 ( 7T 2 - 6ivr - In 8 (In 8 - 4) 



+ q 2 ro 2 ^6ivr(2 + u 2 ) - vr 2 (l + u 2 ) - 24 (u 2 - u + In 2) + In 8 (in 8 + u 2 (4 + In 8)) 
- 2(q 2 - 3t» 2 )(q 2 (l + u 2 ) - 3n> 2 ) {-m + log (1 - u)) 

+ 4 ^q 4 (l + u 2 ) + 9t» 4 (ln 2 - 1) - 3» 2 q 2 (in 2 - 2 + u 2 (ln 2 + 1)) \ In (1 + u) 
+ 3t» 2 (q 2 (l + u 2 ) - 3w 2 ) ln 2 (l + u) 

+ 2 (q 2 (l + n 2 ) - 3» 2 ) (q 2 - 3ro 2 (l + In 2) + 3t» 2 In (1 + u)) In (1 - u) 

'l + u 



+ 6t» 2 (q 2 (l + u 2 ) - 3rx> 2 ) Li 2 



(A.8) 
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